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After a brief introduction to low-x QCD and to resummation approaches, I illustrate the predictions for the gluon Green function and 
splitting function of a recent renormalization-group improved small-x resummation scheme. I argue, on this basis, that the range of 
validity of perturbative calculations is considerably extended in rapidity with respect to leading log expectations. The perturbative high- 
energy exponents are predicted in a phenomenologically interesting range, and significant preasymptotic effects are found. In particular, 
the splitting function shows a shallow dip in the moderate-x region, followed by the expected small-x power increase. 



1 Low-x Physics in QCD 

The description of high-energy hard cross-sections in 
QCD perturbation theory has a long history. Starting from 
the early BFKL prediction of rising cross sections 1 1 1 
at leading logs level (LL) and from the fc-factorization 
framework \ 2\, up to the difficult calculation of the next- 
to-leading (NLL) kernel (3]|4] and to its improvements 
0|i0H|9][To|[n][T2|, a lot of work has been devoted to 
fill the gap between theory and small-x phenomenology. 

On the experimental side, the DIS structure functions at 
HERA 1 13 1 show a rise which is well described by large 
scaling violations arising in low order DGLAP 1 14 1 evolu- 
tion. On the other hand, data for two-scale processes like 
f{Q) -f(Qo) for which the DGLAP picture is inappro- 
priate, seem incompatible with the strong rise predicted by 
LL BFKL evolution. There is thus a need to reconcile the 
two theoretical pictures, which actually differ by the order 
of resummation of large logarithms (log Q 2 vs. log s) in 
the perturbative series. 

On the theoretical side, one expects QCD perturbation the- 
ory at low x to be governed by some effective coupling 
a s (k 2 ) log l/x which risks to be strong for two independent 
reasons. Firstly, because of the logarithmic enhancements 
due to the opening of the smalk gluon phase space. This 
motivates the calculation of next-to-leading corrections to 
the hard-Pomeron intercept, which however turn out to be 
quite large and negative |3,4|, pointing to some instabil- 
ity of the log s hierarchy. Secondly, because energy evo- 
lution - by diffusion or tunneling 1 15 1 - may give rise to 
low values of k, the gluon momentum transfer, thus im- 
plying a drift towards a really non-perturbative asymptotic 
regime ll6H17lfT8lll9l . characterized by soft physics. We 
must control, therefore, higher order subleading terms, and 
analyze the perturbatively accessible rapidity range. 

A tentative answer to the above question marks about the 
perturbative behaviour comes from the (double) resumma- 
tion approaches. The idea is to take into account higher or- 
der subleading terms in both variables (log s and log Q 2 ), 



based on some distinctive features of such terms. For in- 
stance, the "duality" approach |8 9| starts from the ap- 
proximate symmetry between scaling violations and en- 
ergy dependence in order to hint at a resummed form of 
the gluon splitting function. The "renormalization-group 
improved" (RGI) approach 1 5 6 7 1 starts from BFKL evo- 
lution at NLL level, and resums those higher order sub- 
leading terms which are required by the RG behaviour for 
both Q ^> Qo and Qq ^> Q. Other resummation models are 
available also I I I Oil 11121 . The RGI approach is able to pro- 
vide predictions |20| for both the gluon splitting function 
and the gluon density itself, and I will concentrate on that 
approach in the following. 

The BFKL equation is, to start with, an evolution equa- 
tion in the rapidity Y = \og(s / QQq), induced by a kernel 
K(k,k') in k -space, which can be expanded in a s and has 
been calculated up to second order. The LL level is dom- 
inated by (effective) ladder diagrams with gluon exchange 
and the corresponding kernel is scale-invariant and given 
by (a,. =N c a s /%). 



a s Ko(k,k') 



a. 



(k-k') 2 



virtual terms. 



(1) 



The highest eigenvalue of Kq yields the exponent for the 
energy dependence and is co s = 41og2a v , a value which is 
too large compared to HERA data (e. g., 0.55 for a s =0.2, 
compared to an observed value of about 0.2). Similarly, the 
gluon anomalous dimension, given by a s /(0 in » = N — 1 
moment space at one-loop level, acquires higher order sin- 
gularities with positive residues. Roughly speaking, the 
LL approximation grossly overestimates both two-scale 
cross-sections' rise and scaling violations. 

The NLL kernel, extracted in |3]@] after many years of 
work of various groups, yields some qualitatively new fea- 
tures and question marks. Firstly, the corrections to the 
high-energy exponent are negative (which is good), but 
with a quite large coefficient, as follows 



to 



(1-6.47% + ...). 



(2) 
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This points towards an instability of the leading-log s hi- 
erarchy, which is confirmed by the risk of oscillatory be- 
haviour of the gluon density 1 2 1 1 for sizeably different ex- 
ternal scales. Secondly, running coupling effects are taken 
into account, and the scale of the coupling suggested by 
the kernel is (k — k') 2 = q 2 . Therefore, in order to avoid 
divergent integrations, we need to regularize the Landau 
pole, by introducing some cutoff (or freezing parameter) 
k. However, when the energy increases, there is a favoured 
evolution towards smaller values of k = G (k) , by dif- 
fusion or tunneling, and this implies that the asymptotic 
high-energy exponent becomes ^-dependent, i.e., non- 
perturbative (asymptotic Pomeron). Therefore, NLL cal- 
culations raise more questions than they are able to answer, 
and call for the inclusion of higher order terms. 



2 The Renormalization Group Improved 
Approach 

Higher order subleading contributions are not known in de- 
tail. We can argue, however, that a class of them are large, 
and required in order to build single-logarithmic scaling 
violations when either Q or Qq is the leading scale l5l . 
and Q 2 /s or Qq/s is the corresponding Bjorken variable. 
Taking into account such terms is equivalent to resumming 
subleading terms in log s which are known because they 
are leading in log Q 2 and are provided therefore by the 
renormalization group. This idea was used in order to build 
the RGI approach of |7|, where the solution of the homo- 
geneous small-jc equation was studied in detail by the use 
of the ©-expansion method |6| so as to calculate stable 
high-energy exponents and gluon anomalous dimension. 
In recent papers |20| the approach has been extended to 
the full inhomogeneous equation in a slightly different re- 
summation scheme, which is more suitable for numerical 
evaluation. I am thus able to describe results for both the 
two-scale gluon density and for the splitting function in the 
collinear limit. Results for the splitting function have been 
recently obtained in the duality approach also |22|. 

The basic problem considered in 1201 is the calcula- 
tion of the (azimuthally averaged) gluon Green function 
G(Y;k,ko) as a function of the magnitudes of the external 
gluon transverse momenta k = |k|, ko = |ko| and of the ra- 
pidity Y = log^. This is not yet a hard cross section, 
because we need to incorporate the impact factors of the 
probes Il2l l23ll24l . Nevertheless, the Green function ex- 
hibits most of the physical features of the hard process, 
if we think of k 2 , k^ as external (hard) scales. The limits 
k 2 ^> k,Q (£q ^> k 2 ) correspond conventionally to the ordered 
(anti-ordered) collinear limit. By definition, in the ©-space 
conjugate to Y (so that © — dy) we set 



and 

©^(k,ko) =S 2 (k-k ) + y</ 2 k\^4k,k')^4k',k ),(4) 

where J^(k,k') is a kernel to be defined, whose © = 
limit is related to the BFKL F-evolution kernel discussed 
before. 

The precise form of the kernel ,Jf m is given in Ref. H20I 
However, the basic features of the RGI approach is illus- 
trated by a simple observation: in BFKL iteration, all pos- 
sible orderings of transverse momenta are to be included, 
the ordered (anti-ordered) sequence k k\ ■ ■ ■ k n ■ ■ ■ 
ko (k -C k\ ■■■ -C k n •••<?; ko) showing scaling violations 
with Bjorken variable k 2 /s (^/i). Therefore, if only lead- 
ing log A: 2 contributions were to be considered, the ker- 
nel J^ca acting on p- {^j^j would be approximately repre- 
sented by the following eigenvalue function (in the frozen 
coupling limit) 

-^o^a s (—^ + —^ ) ( -+Ai(a>) J +••• 

a s = a s —, (5) 

where Ygg = <%s +Aj(©)) is the one-loop gluon-gluon 
anomalous dimension and we have introduced the variable 
y conjugate to logfc 2 . Note, in fact, that Eq. l|5} reduces to 
the normal DGLAP evolution 1 14 1 in log A: 2 (logA: 2 ,) in the 
two orderings mentioned before, because y+ (1 + y — 
7) is represented by <3 logjt 2 (^log^) at fixed values of x = 
k 2 js (xq = k^/s) in the ordered (anti-ordered) momentum 
region. Note also the ©-dependent shift 1 5 6 7 1 of the y- 
singularities occurring in Eq. (0, which is required by the 
change of scale (k^ versus k 2 ) needed to interchange the 
orderings, i.e., xq versus x. 

We thus understand that the ©-dependence of is essen- 
tial for the resummation of the collinear terms and can be 
used to incorporate the exact LL collinear behaviour, while 
on the other hand, the © — > behaviour of ,J4f m is fixed by 
the BFKL limit up to & (©) terms, so as to incorporate ex- 
act LL and NLL kernels. Such requirements fix the kernel 
up to contributions that are NNL in lnx and NL in ln(2 2 - 
The resulting integral equation to be solved by the defini- 
tion @ is thus a running coupling equation with non linear 
dependence on a s at appropriate scales, and it has a some- 
what involved ©-dependence in the improved kernels J(f w . 
Its solution has been found in |20| by numerical matrix 
evolution methods in k- and x- space |25|, where the typ- 
ical ©-shifted form in the example (0 corresponds to the 
so-called consistency constraint [ 26 27 28 1. Furthermore, 
introducing the integrated gluon density 



Sffl,(k,ko) =[co-Jf a ]-\^o), 



(3) 



xg(x,Q 2 ) = J d\G iso=k2) (logl/x;\klko), (6) 
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A*o=k 

*(6 



(7) 



10 



the resummed splitting function P e s{z,Q 2 ) is defined by 
the evolution equation 



dgj^Q 2 ) fdz n 2sJ x n 2 



(8) 



and has been extracted [20 1 by a numerical deconvolution 
method |29|. We note immediately that P e ff turns out to 
be independent of Icq for Q 2 3> k\, yielding an important 
check of RG factorisation in our approach. 



3 Gluon Green's Function 

Results for G(Y;k,k) are shown in Fig. ^ In addition 
to the solution based on our RGI approach (NLLb) the 
figure also has 'reference' results for LL evolution with 
kernel ^(x^q^K®. The one-loop running coupling (with 
rif = 4), is regularised either by setting it to zero below a 
scale k ('cutoff') or by freezing it below that scale (a s (q 2 < 
k 2 ) = a s (k 2 )). The cutoff regularisation is supposed to be 
more physical since it prevents diffusion to arbitrarily small 
scales and is thus more consistent with confinement - ac- 
cordingly three cutoff regularisations are shown and only 
one frozen regularisation. The k — 0.74 GeV cutoff solu- 
tion is presented together with an uncertainty band associ- 
ated with the variation of x 2 ^ between j and 2, x^ being a 
renormalization-scale testing parameter, such that a s {q 2 ) 
is replaced by a s (x 2 g 2 ) +b a 2 log x 2 in the part of the 
kernel linear in a s . 

Solutions of (0J with an IR-regularised coupling generally 
have two domains I16II17|[T8*1|15I . separated by a criti- 
cal rapidity Y c (k 2 ). For the intermediate high-energy re- 
gion 1 <C Y < Y c (k 2 ), one expects the perturbative 'hard 
Pomeron' behaviour with exponent <o s , 

k 2 G{Y;k,k) ~ -^exp[a) 4 .(a i (^ 2 ))y + A(a,.,y)] , (9) 



and diffusion corrections H30II3 11 132 33 1 parametrised by 
A(a„,Y). Beyond Y c , a regularisation-dependent non- 
perturbative 'Pomeron' regime takes over 



k 2 G{Y;k,k)' 



12 \« 



k 2 



LL : 
NLL B 



(Op 



(10) 



where the non-perturbative exponent (Op satisfies 171 (Op ~ 
(O s (a s (k 2 )) and hence is formally larger than co s (a s (k 2 )). 2 

1 Actually,slightly different values of the scales in G are taken, namely 
G(Y;k + e,k—e) with e = 0. Ik, in order to avoid sensitivity to the dis- 
cretisation of the 5-function initial condition in|4](cf. Ref. 1 20 for a de- 
tailed discussion). 

2 The behaviour 1 1 01 with (Op > 0) s (a s (k 2 )) is a general feature of linear 
evolution equations such as|4] but not of actual high energy cross sections, 
which are additionally subject to non-linear effects and confinement. 
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Figure 1. Green's function calculated with four different in- 
frared regularisations of the coupling, shown for LL and RGI 
NLL ('NLLb') evolution. The bands indicate the sensitivity of the 
k = 0.74 GeV results to a variation ofx^ in the range j to 1. The 
left and right hand plots differ only in their scales. 



The value of Y c depends strongly on k. In the tunnelling 
approximation, it can be roughly estimated by equating 
eqs. (|9j and (II 01 to yield I34II15I . for any given regular- 
isation procedure, 



Yc(k 2 ) 



glog(fr 2 /P) 
(Op — a> s (a s (k 1 )) 



(11) 



(again with 1 + (Op — > 1 for LL), showing an approximately 
linear increase of Y c with log A: 2 . 
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Within this logic, several aspects of Fig.^are worth com- 
menting. The most striking feature of the LL evolution is 
its strong dependence on the non-perturbative regularisa- 
tion, even for rapidities as low as 5. The exact value of 
Y c depends on the regularisation being used, ranging be- 
tween 5 and 10. In contrast, NLLb evolution remains un- 
der perturbative control up to much larger rapidities and the 
NP pomeron behaviour takes over only for Y > 25, where 
the three cutoff solutions start to diverge 3 . Therefore, Y c is 
considerably larger for the resummed evolution, as conse- 
quence of the fact that subleading corrections lower both 
the PT and - even more - NP exponents (see Eq ll li . 

We note that at their respective F c 's the NLLb Green's 
function is an order of magnitude larger than the LL one: 
the subleading corrections increase the overall amount of 
BFKL growth remaining within perturbative control. How- 
ever, large densities G ~ 1 / a s are reached at cosiderably 
larger values of Y, so that saturation effects 1 30 1 are pushed 
towards higher energies. For the reference values k = 4.5 
GeV, a s ~ 0.2 this translates to Y of order 15, close to the 
kinematical limit of LHC, as a very rough estimate. 

In Fig^only a single value of k is considered. The ques- 
tion of NP contributions is summarised more generally in 
Fig. [2] which shows contour plots of the logarithmic spread 
of the four regularisations. Darker regions are less IR sen- 
sitive, and contours for particular values of the spread have 
been added to guide the eye. Here too one clearly sees the 
much larger region (including most of the phenomenolog- 
ically interesting domain) that is accessible perturbatively 
after accounting for subleading corrections. 

So let us now therefore return to Fig. [2 and examine the 
characteristics of the NLLb Green's function in the pertur- 
batively accessible domain, which should be describable 
by an equation of the form l|9}. The first feature to note is 
that the growth starts only from Y ~ 4. This suggests that 
at today's collider energies (implying Y < 6 I35II36I "). it 
will at best be possible to see only the start of any growth. 
This preasymptotic feature is partly due to the slow open- 
ing of small-x phase space 1371 implicit in our a)-shifting 
procedure. 

Once the growth sets in, the issue is to establish the value of 
CO s appearing in ©. This is a conceptually complex ques- 
tion because in contrast to the fixed-coupling case, co s no 
longer corresponds to a Regge singularity. 

There are running-coupling diffusion corrections 
A(a s ,Y) 13013 1I32I33I . whose leading contribution, ~ F 3 , 
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Figure 2. Contour plots showing the sensitivity ofG(Y,k + £,k — 
e) to one's choice of non-perturbative regularisation, as obtained 
by examining the logarithm of the ratio of the regularisations giv- 
ing the largest and smallest result for G. Darker shades indicate 
insensitivity to the NP regularisation, and contours have been 
drawn where the logarithm of the ratio is equal to 0.1, 0.2 and 
0.4. Plot (a) shows the result for LL evolution, while (b) shows 
RGI NLL evolution (NLLb )■ The regularisations considered are 
those of Fig. ^\ 



for this model is 1 20 



3 Renormalisation scale uncertainties of the resummed results are sizeable 
- of the order of several tens of percent for Y > 4 - but seem anyhow quite 
modest compared to the order(s) of magnitude difference with LL. The Xu 
dependence of the LL solution is somewhat smaller than for NLLb- this 
may seem surprising, but at larger Y the LL solution is in the NP domain, 
where non-linearities (in a s ) reduce the dependence. 



A(a s ,Y). 



24 



d log k 2 



(o s (a s (k 2 )) 



Xeffij) 



(12) 



In addition, A(a s ,Y) contains terms with weaker Y depen- 
dences, including Y 2 and Y. Such terms can be disen- 



QCD@Work 2003 - International Workshop on QCD, Conversano, Italy, 14-18 June 2003 



5 



tangled by the method of the ^-expansion |38|. . Since 
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Figure 3. The small-x exponents: the Green's function effec- 
tive exponent (O s is shown to first order in the b-expansion; the 
splitting function exponent (O c is shown together with NP and 
renormalisation scale uncertainty bands, defined in figure^ Also 
shown, for reference, is the result for 0} c using the method of [ 71, 
for b(iif = 4). 



many arguments in the BFKL framework, ll7lll6l[T9"i r71l9l 
29 1, have been given in favour of factorisation, Eq.|H] with 
the small-x splitting function Pg g (z } Q 2 ) being independent 
of the IR regularisation. The most dramatic demonstration 
of factorisation is perhaps in the fact that a numerical ex- 
traction of the splitting function from the Green's function 
by deconvolution gives almost identical splitting functions 
regardless of the regularisation. This is illustrated in Fig.0] 
where the solid line and its inner band represent the result 
of the deconvolution together with the uncertainty result- 
ing from the differences between the three cutoff regular- 
isations. The resulting regularisation dependence is pretty 
small, and at higher Q it diminishes rapidly as an inverse 
power of Q, as expected from a higher twist effect. 
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running-coupling diffusion corrections start only at order 
b 2 , it is possible to give an unambiguous definition of (O s 
up to first order in b, while retaining all orders in a s for non 
running-coupling effects. The result is shown in Fig.[3]as a 
function of a s and as has been found in previous work 1 7 1, 
there is a sizeable decrease with respect to LL expecta- 
tions. 5 Furthermore, the leading diffusion corrections in 
(I12> turn out to be numerically small, about an order of 
magnitude down with respect to the LL result, due to a 
sizeable decrease of the diffusion coefficient over and 
above the decrease already discussed for co s . 



4 Resummed Splitting Function 

The Green's function G(Y;k,ko) has been investigated in 
the collinear limit k ^> ko ~ k also. In such a case the sen- 
sitivity to the IR regularization is much stronger. However, 

4 In the b — > limit, with a s (k 2 ) kept almost fixed, the non-perturbative 
Pomeron is exponentially suppressed I38I9I . so that the fe-expansion can 
also be used as a way of defining a purely perturbative Green's func- 
tion without recourse to any particular infrared regularisation of the cou- 
pling l38l . 

5 A direct comparison with earlier results for £0, [7 J is not possible, be- 
cause they are based on a different definition (the saddle-point of an effec- 
tive characteristic function), which is less directly related to the Green's 
function. Nevertheless, the present results are consistent with previous 
ones to within NNLL uncertainties. 



Figure 4. Small-x NLLb resummed splitting function, com- 
pared to the pure 1-loop DGLAP and the (fixed-coupling) LL 
BFKL splitting functions. The central NLLb result corresponds 
to Xfi = 1, k = 0.74 GeV; the inner band is that obtained by vary- 
ing k between 0.5 and 1 .0 GeV, while the outer band corresponds 
to j <x 2 l <2. 



Several features of the resummed splitting function are 
worth commenting and comparing with previous NLL cal- 
culations, including various types of resummations Q[9] 
1101 . Firstly, at very large-x it approaches the normal 
DGLAP splitting function. The momentum sum rule is 
satisfied to within a few parts in 10 4 . At moderately 
small-x the splitting function is quite strongly suppressed 
with respect to the LL result and shows not a power 
growth, but instead a significant dip (of about 30% rel- 
ative to the LO DGLAP value, for a s = 0.215). Dips 
of various sizes and positions have been observed before 
in 1 10 8 29 1 though this is significantly shallower than that 
found in 1 10 1 at NLL order and similar to that found in the 
©-expansion |20| and in the duality approach |9 22 1. 

At very small-x one finally sees the BFKL growth of the 
splitting function. We recall that the branch cut, present 
for a fixed coupling, gets broken up into a string of poles, 
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with the rightmost pole located at CO c , to the left of the orig- 
inal branch point (co s ), CO s — CO c ~ b 2 ^a^ 3 |7 1. The origin 
of this correction is similar to that of the b 2 ^a^ 3 contribu- 
tions to (Op for cutoff regularisations 1391 . The dependence 
of co c on Q is shown in Fig. [3] together with its scale and 
IR regularisation dependence. It is slightly lower than the 
earlier determination in the ©-expansion (7) 6 . Both deter- 
minations are substantially below co s , as expected. 

To sum up, the resummed gluon density at comparable 
scales k ~ ko stays perturbative for a wide rapidity range - 
including most of the phase space available at next genera- 
tion colliders - due to a suppression of the non-perturbative 
Pomeron and of diffusion corrections (Fig|2j. The ex- 
pected increase with energy starts slowly (FigQ, due to 
small-x phase space effects, and is regulated by the expo- 
nent co s (Fig. |3 which, in the relevant Q-range is roughly 
of size a s (Q 2 ). In the collinear region k^$> ko the gluon 
splitting function can be factored out, shows a shallow dip 
in the moderate small-x range (Fig|4j, and its increase is 
regulated by the exponent co c < CO s , also shown in Fig. [5] 

A realistic prediction of cross-sections requires the inclu- 
sion of impact factors along the lines of |20|, and of the 
quark sector along the lines of [ 19 7|. But we can already 
say that resummed results show interesting preasymptotic 
effects and are, very roughly, closer to low order predic- 
tions than expected. This in turn may provide a prelim- 
inary explanation of the apparent smoothness of small-x 
cross-sections despite the occurrence, in their description, 
of large perturbative coefficients and of various strong- 
coupling phenomena. 
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